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Abstract. Some years ago, one of the authors (MM) revived a concept to which he gave 
the name of single-particle Dirac oscillator, while another (CQ) showed that it corresponds 
to a realization of supersymmetric quantum mechanics. The Dirac oscillator in its one- and 
many-body versions has had a great number of applications. Recently, it included the ana- 
lytic expression for the eigenstates and eigenvalues of a two-particle system with a new type 
of Dirac oscillator interaction of frequency uo. By considering the latter together with its 
partner corresponding to the replacement of uo by — u, we are able to get a supersymmetric 
formulation of the problem and find the superalgebra that explains its degeneracy. 
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1 Introduction 



A Dirac equation with an interaction linear in the coordinates was considered long ago QTJ] 
and revived more recently || with the name of Dirac oscillator as, when reduced to the large 
component, it corresponds to a standard harmonic oscillator with a very strong spin-orbit 
term. 

The concept gave rise to a large number of papers concerned with its different aspects [[J . 
Two of the present authors (MM and CQ) were particularly interested in the symmetry Lie 
algebra 0] that explains the degeneracy of the spectrum of the one-body Dirac oscillator, 
and later one of them (CQ) showed |J that the problem corresponds to a realization of 
supersymmetric quantum mechanics ||. 

One of the authors (MM) and his collaborators were particularly concerned with two- 
and three-body systems with Dirac oscillator interactions, because in the former case it 
gave insight into the quark-antiquark system and the mass formula for mesons [0, |§ f|, 
while in the latter it could be applied to the three-quark system and the mass formula of 
baryons [|10[ PH . 

In all the extensions to more than one particle, the analysis was made in terms of a 



single Poincare invariant equation of the type employed by Barut and his collaborators [12 



The most usual approach to the n-body problem is through the use of n separate equations 



but that satisfy appropriate compatibility conditions |L3L O, 15 



In a recent paper |1(J, Del Sol Mesa and Moshinsky showed that the two approaches 
mentioned in the previous paragraph are equivalent, at least for a two-body problem with 
a modified type of Dirac oscillator interaction of frequency uo. This last problem, with a 
fully Poincare invariant formulation, has, in the centre-of-mass frame, very simple analytic 
expressions for its eigenvalues and eigenstates, and it displays an extraordinary accidental 
degeneracy. In the present paper, we shall prove that when the problem of frequency u is 
taken together with that of its supersymmetric partner, corresponding to the replacement 
of u by — to ( and referred to in the following as the Dirac oscillator with negative frequency), 
a superalgebra explains the degeneracies that appear. 
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Before entering into the objective of the paper, indicated at the end of the last paragraph, 
we note that in the one-body case, Quesne [5[] considered as supersymmetric partners the 
large components of particle and antiparticle wave functions. It is equivalent though 
to consider as supersymmetric partners the large components of the wave functions of 
frequency uo and — u, and it will be this point of view that will be generalized in the present 
paper. 

In section 2, we determine the spectrum and the eigenstates of the two-body Dirac 
oscillator and of its partner with negative frequency. In section 3, we then reformulate 
the two-body problem in the framework of supersymmetric quantum mechanics and obtain 
the superalgebra explaining the spectrum degeneracies. Finally, section 4 contains the 
conclusion. 

2 The two-body problem with a new type of Dirac 
oscillator interaction 

In ref. [[3I|, a two-body problem with a new type of Dirac oscillator interaction was in- 



troduced through a Poincare invariant equation. In the centre-of-mass frame, the latter 
becomes 

'(«! - a 2 ) ■ (p =F ifrA/^^) + Pi + fa] ^ = Ei/j± (2.1) 

where r, p are defined in terms of the coordinates X\, x 2 , and momenta p 1( p 2 of the two 
particles as 

r = x l -x 2 p=\{p l -p 2 ) (2.2) 
while cx.\ and a 2 are the direct product matrices 

— U?)«(S!) *-(S!)«U?) «> 

and similarly for fii, (3 2 , 751, and 752. Note that as we want to simultaneously analyse the 
Dirac oscillators for frequencies u and —u, in eq. (2~T) we have written a =F sign before the 



term proportional to r, and we have appended a corresponding ± superscript to ip. 
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It is now a question of discussing the eigenstates ip of ( |2.1[ ) and the corresponding 
eigenvalues of the energy. We first note that from 
written as 4 x 4 matrices in the form 



, ai, a 2 , /?i/?275i752, /3i + /?2 can be 
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(2.5) 



where for the moment we suppress the ± sign as we first deal with the case of positive 
frequency uj. 



If we now introduce the notation |L6 

a± = (01 ■ P) ± i{oj/2){cr 2 ■ r 
equation ( |2.1| ) can be written as 



{(T 2 -p)±i(u/2)(tT 1 -r) 



(2.6) 
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(2.7) 



The second and third rows of the matrix operator equation ( |2.7| ) allow us to express 
^21, ^12 in terms of if)u, ip22, so that substituting them in the first and fourth rows, we 
obtain a 2 x 2 matrix operator equation for the two components ipn and ip22- Introducing 



then, as in previous work [17], <pi and 02 by the definitions 



1 

V2 



1^22 



{21 
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we finally obtain an equation of the form 



A+ - E 2 2E 
2E B + - E 2 



(2.9) 



where 



A + = (a_ - 6+)(a+ - 6_) = 4cu [S 2 + (S • rj)(S ■ £) + L • 5 

= 4w(5.€)(flf.|j) (2.10) 

5+ = (a_ + 6 + )(a + + 6_) = 4w[jV'-(flf-f7)(5-^)-i-S 

= MS'-tjKS'-O 



(2.11) 



with creation 77 and annihilation £ operators defined by 



77 = ^ ((^/2) 1 /V - i( U /2)-^p) £ = ^ ((^/2) 1 /V + i( U /2)-^p) (2.12) 



and 



N=T)-$ > 



r x p = —i(r] x £) 



while 



!(<7l + 0" 2 ) 



5' 



(<Ti - <7 2 ). 



(2.13) 



(2.14) 



Note that we have added an upper sign + to A, B, 0i, 2 to indicate that they belong 
to the equation where the frequency is u. From ( |2.1|) , we see that the equation for negative 
frequency — uj can be reduced to that for positive frequency u if we replace r by — r, which 
from ( |2~T2D and (fTT3|) implies that 77 -> £ -> -77, iV -> JV + 3, L -> -L. 

Thus for the corresponding problem with frequency — u;, we have the equation 



A" - E 2 2E 
2E B~ - E 2 







where 



A' 
BT 



4co 



S 2 + (S- £){S ■r J )-L-S\= 4u(S ■ rj){S • £) 

N + 3 - (5 ■ g)(S ■ 77) + L ■ S] = 4lu(S' • £)(S" ■ 77). 



(2.15) 

(2.16) 
(2.17) 



It is now a question of finding the states <pf, (f)f and the corresponding energies. This 
is facilitated by the fact that A ± , B^ commute with the total number of quanta N and 
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angular momentum J = L + S. Hence, we could express the states (J)f, (f)f in terms of 
harmonic oscillator states with spin either or 1, i.e., kets of the type 

\N(l, s)jm) = J2 [fa, scr\jm) R Nl (r)Y ¥ (9, ip) X sA (2.18) 

where s — 1, a — 1, 0, —1, or s = 0, a = 0. 

We note that A + and B + commute with one another, as do also A~ and B~ , and 
thus rather than directly look for the eigenvalues of the energy E, we shall first consider 
the eigenstates of A and B and their corresponding eigenvalues, which we shall denote 
by A^ and A^, respectively. 

We shall start our discussion by remarking that for fixed values of N and j, and for 
s = 0, there is only one ket of the form ( 2.18Q , whose parity is (—!)•', namely 



MN) = \N(j,0)jm) (2.19) 

where we introduced the shorthand notation (po{N) for this ket. Note that we indicate the 
explicit dependence of (fo on N, as this value will change in the following section, but we 
do not include a dependence on j as this non-negative integer will remain fixed. 

When the spin is 1, but the parity continues to be (— there is again only one ket, 
which is 

<p' Q (N) = \N(j,l)jm) (2.20) 

where all states of spin 1 will be denoted by <p'. On the other hand, if s = 1 but the parity 
is — (— iy, there are then two kets 

ip' ± (N) = \N(j±l,l)jm). (2.21) 

If A*, A^ can be determined in the various instances discussed in the previous paragraph, 
then the matrix operators in eqs. (|2.9| ), (|2.15|) become purely numerical matrices 

( \±-E 2 2E \ 



2E \f - E 2 



(2.22) 

which will give rise to secular equations for the energy of the form 

(E 2 - Xfj (E 2 - Xf) - AE 2 = 0. (2.23) 
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We shall now proceed to obtain these equations when the states are given by ( |2.19| ), ( |2.20| ) 
and ( |2.21| ), respectively. 

Let us start with the case of positive frequency u and spin s = 0. From (|2.10|) , ( p. 11 



and (|2.19|) , we obtain 



A+ = Xt = AuN 



and thus the square of the energy in ( |2.23j ) can take the values 



E 2 = E 2 = 4 + AuN. 



(2.24) 



(2.25) 



The vanishing value gives rise to the phenomenon we have called a cockroach nest fllSfl . So, 
in the following analysis, we shall disregard it and be only concerned with positive E 2 , such 
as the second expression in ( |2.25|) , which gives an equally spaced spectrum for E 2 . 

Turning now our attention to the state of spin 1 and parity (— 1) J ', we note from Q2.20| ) 
that H 



(N(l,l)jm\L-S\N(l,l)jm) = § + 1) - 1(1 + 1) - 2] (2.26) 
(N(jA)jm\(ri.S)(£-S)\N{j,l)jm) = N + l. (2.27) 



Thus from the middle expression in eqs. ( [2.10| ), ( p. 11] ), we obtain that 



A+ = 4^(iV + 2) A+ = 



and the relevant square of the energy takes then the form 



E 2 = A + Au{N + 2). 



(2.28) 



(2.29) 



For the case of spin 1 and parity — (— 1) J , the kets are \N(j ± 1, l)jm), and the non- 
vanishing matrix elements of L ■ S are given by ( [2.26|) , while those of (rj • S){£ • S) were 



obtained in ref. 191. Because of the existence of two kets instead of one as before, the 



operators A + and B + now become 2x2 numerical matrices. Their eigenvalues 



A+ = A+ = 4^(iV + 2) 



(2.30) 
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correspond to the eigenstate 

, * m ( (j + l)(N + j + 2) \ 1/2 / j(N-j + l) \ 1/2 

^ l(iV)= ( (2j + !)(* + 2) J V9+(jV) + ((2 J + l)(iV + 2)J ^ (2 ' 31) 

while the eigenvalues 

A+ = 4cu(iV + 2) A^ = (2.32) 

correspond to the eigenstate 

^ (JV)= ( (2j + l)(i V + 2) J MJ0-( (2j + 1)^ + 2) J (2 ' 33) 
In both cases, the relevant square of the energy is 

E 2 = 4 + 4u(N + 2) (2.34) 



so from ( 2.29| ) and ( |2.34| ), we see that we have the same E for all s = 1 states. 



In the case where the frequency is — u, the analysis is entirely similar, and only the 
results will be given. For s = and parity (— l)- 7 , we have 

\~ = Q X^=Auj(N + 3) E 2 = 4 + 4u(N + 3). (2.35) 

For s = 1 and parity (— iy , we get 

= 4a;(iV + 1) A^ = E 2 = 4 + Auj(N + 1). (2.36) 

Finally, for s = 1 and parity — (— 1) J , we have two possibilities 

A" = A^ = 4c;(iV+l) E 2 = 4 + Auj(N + 1) (2.37) 
A~ = Aco(N + l) A,7 = E 2 = 4 + Auj(N + 1) (2.38) 

so again, for all situations when s = 1, we get the same E 2 . 

The eigenkets for frequency — uj will be denoted with a bar above, i.e., as 7p (N), yf (N), 
Tp'^N), Tp' 2 (N). An analysis similar to the one carried out between eqs. (|2.19|) and Q2.33j ) 
shows that 

Tp (N) = <p (N) Jp' (N) = ip' (N) (2.39) 
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while (pi(N), f' 2 {N) are given by 

/(j + l)(iV-j + l)\ 1/2 



(2j + l)(N + l) 
j(N + j + 2) ^ 



,'an)J ^ n :j: 2 1 ) 1/2 ^n) (2.40) 



1/2 



(2j + l)(N + l) i 

*(j + i)(N-j + iy 



1/2 



<pf_(N) (2.41) 



V(2j + i)(jv+i); " +v 7 v (2j + i)(iv + i) , 

with (^o(iV), ^o(iV), (f/±(N) given by (pl9|) , (^20]) , and ( ^211) , respectively 

Note that whenever the eigenvalue of A is different from zero, that of B ± vanishes, 
and viceversa. From the above analysis, it therefore follows that (4a;) _1 (A + + B + ) has 
eigenvalue N if s = 0, and N + 2 if s = 1, while for (4a;) -1 (A~ + 5 _ ), it is respectively 
N + 1 for s = 1, and N + 3 for s = 0. All these eigenvalues are displayed in table |TJ, where 
besides we give the spin and parity of the corresponding eigenkets. 

The combined spectrum of the square of the energy, in units 4a;, for both frequencies u 
and —uj is illustrated in figure [l], where below the level columns we indicate the value of 
the spin and the sign of the frequency. The spectrum shows an extraordinary degeneracy, 
which we intend to explain in the next section through consideration of supersymmetry 
that eventually leads to an appropriate superalgebra. 

3 Supersymmetry and superalgebra for the two-body 
system with a new type of Dirac oscillator interac- 
tion 

As we showed in the previous section that all states of spin 1, regardless of whether they 
have parity (— iy or — (— l)- 7 , give rise to the same eigenvalue N + 2 for A + + B + , and N +1 
for A~ + B~ , it is convenient to group them as vectors in the notation 



( v'o(N) \ 
cp[(N) 

V <AW J 

( %(N) \ 
Tp[(N) 

{ &{N) ) 



(3.1) 
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where the components in the unbarred and barred vectors were all given explicitly in sec- 
tion 2. 

We remarked previously that (4a;) -1 (A + + B + ), when applied to (p (N), ip'(N), gives 
the eigenvalues N, N + 2, respectively, while (4a;)- 1 (A~ + B~), applied to Tp Q (N), Tp'(N), 
gives N + 3, N + 1. In analogy to the superstate discussed in ref. || for the one-particle 
problem, let us now introduce a superstate for the two-body problem, defined by the ex- 
pression 



(3.2) 



( MN) \ 
T?(N-1) 
cp'(N-2) 
\lp (N-3) J 

with (p , Tp , cp', Tp' defined by ( p. 19 ), ( 2.39| ), and (|3.1| ), respectively. From the discussion of 
the previous section, and the definitions (|2.10| ), ( |2.11|) , (|2.16| ), and (|2.17|) of A ± , B ± , it is 
clear that \l/ is an eigenstate with eigenvalue N of the following supersymmetric Hamiltonian 



( H+ \ / <p (N) \ 

H Ip'(N-l) 

H + (p'(N-2) 

V H- ) V^ (iV-3) J 



N 



( MN) \ 
cp'(N-l) 
cp>(N - 2) 
\<p (N-3) J 



(3.3) 



where 



H~ 



(Au)- 1 [A + + B+) = (S -£)(S ■ 77) + (S' ■ 77) (S' ■ £) 
(4c)" 1 (A- + B~) = (S ■ r,)(S • £) + (S' • £)(S' • 77). 



(3.4) 
(3.5) 



The problem now resides in finding the corresponding supercharges, i.e., the counter- 
parts of the operators Q, of the one-body case given in ref. f|. We shall retain the 
notation Q, Q\ with the latter being the Hermitian conjugate of the former. We would 
like of course to recreate the relation (3.2) in ref. ||, i.e., 



{Q,Qt} =QQ^ + Q^Q = H 



(3.6) 



but also that Q, when applied to would relate <po(N) with (p'(N — 1), ip'(N — 1) 
with (f'(N — 2), and the latter with Tp (N — 3). 
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It is obvious that Q must contain the annihilation operator but in a form that should 
commute with the total angular momentum J = L + S,asj,m are the same for all states 
in This last point can only be achieved if we consider the scalar products of £ with 
the spin operators of the two particles or, better still, with their sum and difference S, S', 
defined in (p. 14 ). Thus Q must be a 4 x 4 matrix of the type appearing in (|3.3|) for H, but 
whose elements depend on (<S-£), (S' •£). We note though from standard Racah algebra f20| 



that S' can connect only states of spin with those of spin 1 or viceversa, while S can only 
connect states with the same spin. As ipo, Tp have spin while <p', Tp' have spin 1, this 
immediately suggests that 



Q = At ■ £ A f = 
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From the relation 



we immediately get 



[S ■ S)(S> ■ = (S ■ rj)(S' ■ ry) 



{Q,Q} = {Q t ! Qt} = {g,gt} = jfJ 



from which it follows that 



[H,Q]= H,Q 



0. 



(3.9) 
(3.10) 
(3.11) 



We have therefore found for the two-body problem the su(l/l) superalgebra 0] charac- 
teristic of supersymmetric quantum mechanics , and analogous to the one-body relations 
obtained in eq. (3.2) of ref . @. 

As a last point, we would like to obtain the superalgebra responsible for the degeneracy 
we observe in figure 1 for our two-body problem with a Dirac oscillator interaction. For 
this purpose, we note that A, At are vector matrices of components 



A,, A] 
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1,2,3 



(3.12) 



and in the appendix we prove that the latter satisfy the following anticommutation relations 

{A,, Aj} = {At, A]} = {\, A]} = Sn I (3.13) 
where I is a 4 x 4 unit matrix. Equation ( |3.13 ) means that Aj, Aj- behave as fermion 



creation and annihilation operators. 

On the other hand, rji, £j are boson creation and annihilation operators, which satisfy 
the following commutation relations 



[VhVj] = 







{ihvi\ = 5 



'J 



(3.14) 



and of course, from their very definition, the boson r/j, £j and fermion At, Aj operators 
commute among themselves. 

The supersymmetric Hamiltonian, in the form ( |3.10 ), can now be written as 



H = ( Vt A t ) (&Aj) + Aj) ( Vi Ai) 



(3.15) 



where repeated indices are summed over i, j = 1, 2, 3. By using all the above-mentioned 
commutation and anticommutation relations, and the property 
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it is straightforward to show that 

( ri t + S 2 

H = rj-^l+ A 1 "- A = 



V 



o 






• £ + 3 - S 2 



















r/^ + 3-5 2 J 



(3.17) 



It is then obvious that all the operators 
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Uij = At£j 



(3.18) 
(3.19) 



commute with H. The operators (|3.18|) and (|3.19| ) are the even and odd generators of a 
symmetry superalgebra respectively, as they satisfy the following commutation and anti- 
commutation relations 



[Cij , C 



(3.20) 
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[Cij, Ci'ji] 


= Sji'Cij* Sij'Ci'j 




(3.21) 


[Cij, Ti'ji] 


fiji'Tiji [Cij, Ui'ji~\ 


—5ij/Ui'j 


(3.22) 


[Cij, Ti'ji] 


= —5ij'Ti'j [Cij, U^ji] 


= 5ji'Uij> 


(3.23) 




= {U ij ,U i , j ,} = 




(3.24) 


{T^, Ui>j/} 


= dji'Ciji + dij'Ci'j. 




(3.25) 



The latter are the defining relations of a u(3/3) superalgebra ||2T | . 

Such a superalgebra fully explains the accidental degeneracy of the levels of the two-body 
system. It indeed contains some operators giving rise to transitions between degenerate 
eigenstates of a given component of the supersymmetric Hamiltonian ( |3.17| ) (namely, the 
even generators Cij and Cij), as well as some operators generating transitions between 
degenerate eigenstates of different components of the same (namely, the odd generators 
and Uij). 

One should remark that the su(l/l) superalgebra of supersymmetric quantum mechan- 
ics is embedded into u(3/3), since the operators H, Q, and Q\ defined in ( |3.17|) , ( |3.7p , 
and (|3.8| ), can be rewritten as 



h = E( c « + cy Q = Y,u« Q ] = J2 T a- ( 3 - 26 ) 

i i i 

4 Conclusion 

The objective mentioned at the beginning of the paper of finding the symmetry superalgebra 
of the two-body system with a new type of Dirac oscillator interaction has been achieved. 
Moreover, we did identify such a superalgebra with u(3/3), and we did show that it contains 
the su(l/l) superalgebra of supersymmetric quantum mechanics as a subsuperalgebra. 

This u(3/3) superalgebra plays the same role for the present problem as the algebra u(3) 
for the standard harmonic oscillator. For the latter, it is well known that by adding the 
operators 77^, i, j = 1, 2, 3, to the w(3) generators rji^j, i, j = 1, 2, 3, one obtains a 
dynamical Lie algebra of the type sp(6,R) p2| . 

Can one find in a similar way a dynamical Lie superalgebra for the two-body system 
with a new type of Dirac oscillator interaction? The answer is of course yes, as one only 

14 



has to add to the u(3/3) generators Cij, Cy, Tij, Uij, i, j — 1, 2, 3, of eqs. ( |3.18| ), ( |3.19| ), 
the operators mentioned in the previous paragraph, plus operators of the form ^A]-, ^Aj, 



AJAJ, AjAj-, i, j = 1, 2, 3. It is indeed straightforward to show that the whole set of 
operators generate a dynamical superalgebra of the type osp(6/6,M) 
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Appendix. Some relations for spin operators 



By using the well-known relation 

(TiCTj = 5ij + ie ijk o-k (A.l) 

for the Pauli spin matrices, it is easy to check that the spin operators S and S', defined 
in (|2.14j) , satisfy the following relations 

S^Sj -\- SjS^ 



S^Sj ~\~ SjSi 



SjSi 








6ij. 



(A.2) 
(A.3) 
(A.4) 



Then eqs. (|3.9| ) and ( |3.16|) directly follow from ( |A.2|) and (|A.4j ), respectively. 

The proof of eq. ( |3 . 1 3|) is also straightforward. Considering first the second anticommu- 
tator in (|3.13| ), we obtain 
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by successively using (|A^) , and ( [PP . Similarly, we get {A^A^} = 

first anticommutator in (|3.13|) . Finally, the last anticommutator follows from 
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(A.6) 



S^Sa + SjSi J 
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Tables and table captions 



Table 1: The eigenvalues of the operators A + , B + , A~, B~ , and of their sums A + + B + , 
A~ + B~ , are indicated in units of 4w. The spin s = or 1 and parity V = (— iy or — (— iy 
of the states are also shown in the first two colums. 



s 


V 




A + 


B + 


A~ 


B~ 


A + + B+ 


A- + B- 





(- 


-iy 





N 





N + 3 


N 


N + 3 


1 


(- 


-iy 


N + 2 





N+l 





N + 2 


N+l 


1 


-(- 


■iy 





N + 2 





N + 1 


N + 2 


N+l 


1 


-(- 


iy 


N + 2 





N + 1 





N + 2 


N+l 
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Figure captions 



Figure 1: Square of the energy levels in units Auo for the two-body Dirac oscillator for 
frequencies +u and —ou, with the sign indicated in the abscissa as also the spin values s = 
or 1. The action of the supercharge operators Q, is also shown. 
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